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Abstract 

In the hot, dense plasma of solar and stellar interiors, 
the Coulomb interaction is screened by the surrounding 
plasma. Although the standard Salpeter approxima- 
tion for static screening is widely accepted and used in 
stellar modeling, the question of dynamic screening has 
been revisited. In particular, Shaviv and Shaviv apply 
the techniques of molecular dynamics to the conditions 
in the solar core in order to numerically determine the 
dynamic screening effect. By directly calculating the 
motion of ions and electrons due to Coulomb interac- 
tions, they compute the effect of screening without the 
mean-field assumption inherent in the Salpeter approx- 
imation. Here we reproduce their numerical analysis 
of the screening energy in the plasma of the solar core 
and conclude that the effects of dynamic screening are 
relevant and should be included in the treatment of the 
plasma, especially in the computation of stellar nuclear 
reaction rates. 

Keywords equation of state, nuclear reactions, nucle- 
osynthesis, abundances, plasmas. Sun: general 



1 INTRODUCTION 

Screening of nuclear reactions in stellar interiors is a 
hotly debated i ssue. The standard Salpeter treatment 
(|Salpeter II1954I) was a necessary improvement in early 



Katie Mussack 

Institute of Astronomy, University of Cambridge, Cambridge, 
CBS OHA, UK 

Los Alamos National Laboratory, X-2, MS PS65, Los Alamos, 
New Mexico, 87545, USA 

Werner Dappen 

Department of Physics and Astronomy, University of Southern 
California, Los Angeles, California, 90089, USA 



equation-of-state work. Many faithful fans consider this 
as evidence that the traditional screening formulation 
must be used for solar models. However, the incom- 
patibility of models generated with the recently revised 
solar abundances and helioseismic results highlights the 
need to re-examin e the physics used to d evelop and ana- 
lyze solar models (jAsplund et al. 1120091 ). In this paper, 
we re-examine screening in the solar core and present 
evidence that dynamic effects must be considered when 
examining screening in stars. 

I.l Electrostatic Screening 

Under the extreme temperatures and densities of the 
solar core, the plasma is fully ionized. The free electrons 
and ions interact with a Coulomb potential energy 
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In this Coulomb system, nearby plasma is polarized by 
each ion. When two ions approach with the possibility 
of engaging in a nuclear reaction, each ion is surrounded 
by a screening cloud. Each ion is attracted to the elec- 
trons and repelled by the protons in its partner's cloud. 
The combined effect of the particles in the screening 
clouds on the potential energy of the pair of ions is 
referred to as screening. This electrostatic screening 
effect reduces the standard Coulomb potential between 
approaching ions to a screened potential which includes 
the contribution to the potential from the surrounding 
plasma. The reduced potential enables the ions to tun- 
nel through the potential barrier more easily, thus en- 
hancing fusio n rates. As illustrated by other authors at 
this nieeting ( Dappen 20091 Baturin 2009 : Straniero 



this nieetmg (lUappen ll^UUijt Ihsaturm M/UUa Ibtraniero 
2009t lYusof 1 120091 ) ■ understanding screening is impor 



ta nt in solar and s tellar modeling 



Salpeter I (jl954( ) derived an expression for the en- 



hancement of nuclear reaction rates due to electron 
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screening. By solving the Poisson-Boltzniann equation 
for electrons and ions in a plasma under the condition 
of weak screening ((/)intoiaction << ksT), Salpeter ar- 
rived at an expression for the screening energy that is 
equivalent to that of the Debye-Hiickel theory o f dilute 
solutions of electrolytes (jPebve fc ffiickellll923h . 
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where the Debye length, Ad, is the characteristic 
screening length of a plasma at temperature T with 
number density n. 
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1.2 Dynamic Screening 



Although Salpeter's approximation for screening is 
widely accept ed, several pape r s ove r the last few 



decades Ce.g. IShayiv fc Shaviv I Il996l : ICarraro et al 



IQSSt Weiss et al. 2001 ) have questioned either the 
derivation itself or the validity of applying the ap- 
proximation to hot, dense, Coulomb systems like the 
plasma of the solar core. Various work d eriving alterna- 



tive f o rmulae for electrostatic s creening ( Carraro et al. 
1988t lOpher fc Qpher I [iooot IShaviv fc Shaviv I Il996 



Shavchenko Il999[ Lavagno fc Quarati l200d Tsvtovich ^ 
200C ) were refuted in subsequent papers (see|Bahcall et al. 
20021 for a summary of arguments in Salpeter's de- 
fense). However, the question of dynamic screening 
remains open. Dynamic screening arises because the 
protons in a plasma are much slower than the elec- 
trons. They are therefore not able to rearrange them- 
selves as quickly around individual faster moving ions. 
Since nuclear reactions require energies several times 
the average thermal energy, the ions that are able to 
engage in nuclear reactions in the Sun are such faster 
moving ions, which therefore may not be accompanied 
by their full screening cloud. Salpeter uses the mean- 
field approach in which the many-body interactions are 
reduced to an average interaction that simplifies calcu- 
lations. This technique is quite useful in calculations 
that rely on the average behavior of the plasma. How- 
ever, dynamic effects for the fast-moving, interacting 
ions lead to a screened potential that deviates from the 
average value. The nuclear reaction rates will there- 
fore differ from those computed with the mean-field 

ap proximation. 

Shaviv fc Shaviv I (Il996l) used the method of molecu- 



does it assume a long-time average potential for the 
scattering of any two charges, which is necessary in the 
statistical way to solve Poisson's equation to obtain the 
mean potential in a plasma. Shaviv and Shaviv at- 
tribute the differences between their simulations and 
Salpeter's theory to dynamic effects. Since their claims 
have been met with skepticism, we have conducted in- 
dependent molecular-dynamics simulations to confirm 
the existence of dynamic effects. 

The viewpoint presented by Bahcall et al., 2002 can 
be summarized with their statement "There is only one 
right answer, but there are many wrong answers." Al- 
though we agree that equation ([2]) is the right answer to 
the question of static screening, we contend that this is 
not the right question to ask. All arguments in favor of 
Salpeter's formulation rely on a mean-field treatment, 
an assumption that must be tested before it is imple- 
mented. The work presented in this paper addresses 
the more appropriate question "is the mean-field ap- 
proach applicable in stellar plasma?" Our work will 
show that there are deviations from the mean field in 
the case of p-p reactions in the solar core and that dy- 
namic screening should be considered in order to obtain 
a more accurate representation. 



2 METHOD 

llow can we test the qu estion of mean-field theor y's rel- 
evance in solar plasma? Shaviv fc Shaviv ( 19961 ) devel- 
oped a method that relies on the techniques of molec- 
ular dynamics to model the behavior of solar plasma 
without using mean-field assumptions. Their simu- 
lations show deviations from mean-field theory that 
would lead to changes in nuclear reaction rate calcu- 
lations. Their claims have been met with skepticism, 
so we replicated and analyzed their work in order to 
resolve the issue. In our previous wo rk (|Mao et al. 



lar dynamics to model the motion of charges in a plasma 
under solar conditions in order to investigate dynamic 
screening. The advantage of the molecular-dynamics 
method is that it does not assume a mean field. Nor 



,2004, 2009; Muss ack etallbood lioOTl) . we examined 
the methods and assumptions used in Shaviv and Sha- 
viv's work, including their treatment of the long-range 
Coulomb force, the effective quantum potentials, and 
the system size. We did not find any problems with 
their techniques. Furthermore, we confirmed that the 
mean-field theory does not adequately describe the be- 
havior of the plasma. 

Here we show that the screening energy of two inter- 
acting protons in our simulation depends on the rela- 
tive kinetic energy of the pair. We also determine the 
dynamically screened interaction potential energy and 
discuss the significance of this result. 

In order to numerically determine the effect of dy- 
namic screening on p-p reactions, we modeled a 3- 
dimensional box of 500 protons and 500 electrons 
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interacting via the Coulomb potential. The efFec- 
tive pai r pote n tials derived for a hydrogen plas r na by 



Barker I (|l97lh : lDeutsch I (|1977[ ): iDeutsch et al. I (|1978 



19791) were employed to include quantum corrections 



The temperature and density of the solar core (T = 
1.6x10^ K, p = 1.6x10^ kg/m^) were used to deter- 
mine the velocities and density of the particles in the 
box. A thermostat was implemented to maintain con- 
stant temperature. Periodic boundary conditions and 
the minimum-image convention were applied. The ve- 
locity verlet alg orithm followed th e time evolution of 
the system. See Mao et al. ( 20091 ) for more details. 

The screening energy was calculated for pairs of ap- 
proaching protons in the following way. For each pro- 
ton, we designated the nearest approaching proton as 
its partner. Then we tracked each pair of protons 
through their approach and subsequent retreat. At 
the point of closest approach, we recorded the separa- 
tion, Vc, and the kinetic energy of the pair, Ekinetic{fc)- 
When the pair was separated by a sufficiently great dis- 
tance, we recorded the kinetic energy of the pair, 
Ekinetic{Rf), and stopped tracking that pair (for our 
simulations, Rf = 2 as where is the Bohr radius). 
At this point, we designated a new partner and repeated 
the tracking process. 

The screening energy of each pair was computed 
from the difference in energy at Tc and at Rf 



(4) 



This represents the energy exchanged between a pair 
and the surrounding plasma. Equation 2] can be ex- 
panded as 



£'kinGtic('^c) H ) - ( i?kinctic(i?/) 



Ri 



(5) 



3 RESULTS 



A key ingredient for the screening energy in equation 
[5] is the difference in the kinetic energy of a pair when 
partner protons are far apart and when they are at their 
closest separation. 



Ai?kinotic — £^kinctic(-R/) — -^kinetic (''c) ■ 
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In figure [U we show the average change in kinetic en- 
ergy of approaching pairs for each distance of clos- 
est approach. For comparison, we also plot the bare 
Coulomb potential and the statically screened Coulomb 
potential as a function of separation. We see that at 
the distance of closest approach, the average kinetic 



energy exchanged between a pair of protons and the 
plasma closely matches the statically screened poten- 
tial. This confirms that Salpeter's static screening can 
successfully describe average properties of the system. 
However, we can see the dynamic effect on screen- 
ing when we sort the pairs of particles by relative ki- 
netic energy. Figure [5] shows the average energy gained 
from the plasma by pairs of protons with a given far- 
apart kinetic energy. This is quite different from the 
Debye-Hiickel screening energy calculated for the aver- 
age closest-approach distance of pairs in each kinetic 
energy bin. We see that pairs of protons with a ki- 
netic energy less than the thermal energy of the plasma 
gain more energy from the surrounding plasma than 
the mean-field result, while pairs with a kinetic energy 
greater than the thermal energy gain less energy than 
the mean-field result and even tend to lose energy to 
the plasma. From this plot, we can estimate the screen- 
ing energy of the p-p reaction at the Gamow energy of 
4.8 kT. 

In order to quantify the dynamic effect on screening 
in the plasma of the solar core, we have split the total 
interaction energy into the Coulomb contribution and 
the interaction energy from the plasma. 
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For the mean-field treatment. 



E. 



screen, mf 
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(l_exp(-r/Az,)) 
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Because there is no formalism to compute the dy- 
namic effect analytically, we use the simulation results 
to determine -Escroen for dynamic screening. One key 
difference from the static screening expression is that 
the dynamic screening energy is a function of both pair 
separation and relative velocity. As before, we split the 
total interaction energy into the Coulomb and screening 
cloud contributions. 



U{r,v) = ?7coulomb(r) - £^scrccn,dyn 



(9) 



For comparison, the dynamic screening energy at the 
distance of closest approach can be written in a form 
similar to the static screening energy 

£'scrocn,dyn(?'c, w) = — (1 - exp{-rc / Ad{v))) (10) 

rc 

by including a new velocity dependent Debye-like ra- 
dius, Ad{v). Figure [3] shows the form of A£){v) deter- 
mined from the simulations. 
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4 DISCUSSION 
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Fig. 1 Dependence of the average kinetic energy change 
<Ai5kinetic> on the closest distance rc, compared with the 
Coulomb potential and screened Coulomb potential. 



This work confirms that screening in the hot, dense 
plasma of stellar cores depends on the relative velocities 
of the interacting ions. The Debye-Hiickel screening 
energy is only valid for describing average properties of 
the plasma. Since faster ions are more likely to engage 
in nuclear reactions than thermal ions, the mean-field 
treatment does not provide an accurate representation 
of this velocity-skewed phenomenon. In fact, the fast 
pairs of ions tend to lose energy to the plasma instead 
of gaining energy from it which would reduce nuclear 
reaction rates instead of enhancing them. Solar and 
stellar models should be adjusted to account for this 
dynamic screening effect. 

Currently, there is no formalism to compute dynamic 
screening analytically. This paper is intended to pro- 
vide insight into the difference between the Salpeter 
formafism and the numerically determined dynamic 
screening. A detailed calculation of the dynamic screen- 
ing correction to the p-p reaction rate in the solar core 
based on this numerical work is underway for a future 
publication. However, these numerical calculations will 
not minimize the need for an analytical formalism for 
dynamic screening in order to generalize the results to 
other temperatures, densities, compositions, and reac- 
tions. Only then can dynamic screening be encorpo- 
rated consistently in solar and stellar models. 
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Fig. 2 Dependence of the screening energy from the sim- 
ulation (squares) on the far-apart kinetic energy El^^^^^^. 
The Debye-Hiickel screening energy computed at the 
averege closest-approach distance of pairs of protons with 
a given far-apart kinetic energy is shown (circles) for com- 
parison. 
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Fig. 3 Modified Debye length for dynamic screening at the 
distance of closest approach obtained from simulations. 
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